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(57) The present invention relates to radio commu* 
nication antennas and in particular relates to multi-user 
detection for CDMA antenna array receivers. A reverse 
link (mobile to base) of a cellular direct sequence code 
division multiple access system is provided which in- 
creases system capacity whilst avoiding system com- 
plexity. The invention provides a method of signal de- 
tection which is based on the partitbning of users into 



equivalence classes by beam-forming with subsequent 
standard multi-user detection techniques applied within 
each class. Computation techniques are simplified and 
processing time can be considerably reduced. Further- 
more, the capacity of a conventional receiver may be 
improved by employing multi-user detection techniques. 
These techniques are particularly useful for systems 
which involve large numbers of users. A method of op- 
eration is also disclosed. 
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Description 

FiELD OF THE INVENTiON 

5 [0001] The present invention relates to radio cellular mobile communications and in particular relates to multi-user 
detection for code division multiple access antenna array receivers. 

BACKGROUND OF THE INVENTION 

10 [0002] To cope with the increasing demand for cellular mobile communications it Is necessary to find ways to increase 
system capacity on the reverse \\r\k whilst avoiding system complexity. 

[0003] The cellular mobile communications IS-95 standard describes the use of direct sequence code division mul- 
tiple access (CDMA) techniques. In such systems, each user is allocated a distinct pseudo-noise (PN) code. The signal 
from each user is multiplied by a respective code before transmission to the base station. All users transmit using the 
IS same radio frequency carrier. The signals from different users will arrive asynchronously due to their different locations 
and signals from each user may arrive asynchronously due to multipath propagation. 

[0004] Figure 1 is an algebraic representatbn of a CDMA communications link. The vector ^contains A/ consecutive 
binary data symbols for P users. When these symbols are transmitted, they are subject to multipath distortion. This 
causes the receiver to observe J versions of each transmitted symbol, which arrive at different times. This effect is 

20 defined mathematically by two matrices. Multiplying d by the matrix T repeats each symbol J times. The size NPJ 
matrix A is diagonal. Its diagonal elements are the positive square roots of the received multipath fading powers for 
the /VP signals received on J paths. This results In the received signal being characterised as the product ATd The 
size A/PJx NPJ matrix (RA)D represents the combined effects of beamforming and of pseudo-noise coding and de- 
coding, where L is the CDMA processing gain. The size A/PJxA/FJ matrix R Is Hermitian [Hom92. pi 69] (R.A. Horn 

25 and C.R. Johnson, "Matrix Analysis", Cambridge University Press, Cambridge (UK), 1 992.) The size NPJxNPJ matrix 
D is diagonal [Hom92. p23]. The quantity y = (R/L)DATd + z/L represents the processed received signal plus back- 
ground noise (the size NPJ vector z/L). 

[0005] The vector y may be subject to conventional bit detection techniques, e.g. as described in [Proakis95] (J.G. 
Proakis, " Digital Communfeattons (3rd Ed)", McGraw-Hill. 1995). Alternatively, the capacity of such a system may be 
30 improved by employing multi-user detection (MUD) techniques in which information about multiple users is used to 
detect a desired user. Another way of increasing the capacity of the system is by employing a steerable beam antenna 
array at the base station. This enables the multiple access interference (MAI) between users transmitting from distinctly 
different bearings to be reduced. However, the MAI between users transmitting from a similar bearing may not be 
reduced. 

3S [0006] It is known that by using multi-user detection [Moshavi 96] (S.Moshavi, "Multi-user Detection for DS-CDMA 
Communrcatlons", IEEE Personal Comms Mag, Vol 34(10), October 1996, pp1 24-35) or antenna array receivers 
[Naguib94] (AF. Naguib, A. Paulraj and T. Kallath, 'Capacity Improvement with Base Station Antenna Arrays in Cellular 
CDMA", IEEE Trans Veh Tech. Vol 43(3), August 1 994, pp 691 -7). bit error rates considerably lower than those provided 
by the conventional detector can be achieved for the reverse link of a cellular direct sequence code division multiple 

40 access system. Results displayed in Figures 8a,b-1 6a,b confirm these findings for both Additive White Gaussian Noise 
(AWGN) channels and Rayieigh fading single and multipath channels. Initial beamforming followed by multi-user de- 
tectbn can further decrease bit error rate for these channels (see Figures 8b-16b), and hence increase capacity, but 
at a cost in complexity. Figure 3 shows a generic receiver employing beamforming foltowed by multi-user detection. 
The system is necessarily complicated and computational requirements are high. 

4S [0007] Several detection systems exist to provide an estimate of a vector of transmitted bits d, given an output y. 
The detectors under consideration are the conventional (single user) detector and four multi-user detectors, namely 
the linear decorrelator detector, the linear minimum mean square error (MMSE) detector, the non-linear decision feed- 
back decorrelator and a form of subtractive interference cancellation (also non-linear). The latter four are usually de- 
scribed for the case of single path transmission (J=1 and T=l, where / is the identity matrix [Horn92.p6]) and in the 

so absence of fading (D=r/), that is. when j[=(R/L)Ad + z/L 

[0008] In the simplest system, single-user detection is employed. The signal from a particular user is detected by 
correlating the received signal, which is a sum of signals from all transmitting users, with the PN code of the user. The 
matched filter detector estimates the transmitted bits according to the signs of the real parts of the received output jf. 
The signals from other users interfere with the desired signal and the system capacity is limited by multiple access 

ss interference. 

[0009] The linear decorrelator detector employs ah inversion of an estimate of the matrix R/L in order to estimate 
Ad. where L is the processing gain. A positive definite Hermitian estimate of /?, say R\ is cateulated using knowledge 
of the PN codes, delays, phases and array signatures of the P users. Then R'A is inverted using the Cholesky decom- 
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position [Horn92, p407] and the transmitted bits are estimated according to the signs of the real parts of the components 
of (R'A)'^^, Estimates of the received signal powers are not required. 

[001 0] The MMSE detector Is represented by an /VP x NP matrix C which minimises the following expression: 

[0011] A closed fomn expression for Ccan be determined following a method suggested by [Honig95] ( M. Honig, U. 
10 Madhow and S. Verdu. "Blind Adaptive Multiuser Detection'. IEEE Trans. Info. Theory, Vol. 41(4), July 1995. 
pp954-960). The result for C is obtained as: 



,5 C^E=((fVL)^(o^/Li)A 

where denotes the background noise variance and A'^ denotes the matrix inverse operation [Hom92, pi 4]. The 
MMSE detector takes into account both the background noise and the received signal powers. In general, the MMSE 
detector does not enhance the noise as much as the decorrelator and so provides a better bit error rate. Estimates of 
20 the received signal powers and the level of background noise are required. 

[0012] The decision feedback decorreiator makes bit decisions in the order of decreasing received signal powers. 
Hence these powers need to be estimated. It employs a Cholesky decomposition to factor the positive definite Hermitian 
matrix R into F^F, where FIs a lower triangular matrix and is the Hermitian adjoint or transpose [Hom92,p6] of F. 
The filter JL (F^*^ is applied to the sampled output y to yield: 

25- 

VL {F/jL)Ad^jL(F^)'' zJL 

[0013] In practice, B can be estimated and hence F. As F/-/L Is tower triangular, the /r-th component of 7l (/=^*^^ 
so does not contain a multiple access interference term for any other bit /c'> fc So the 0-th component does not contain 
an MAI term due to any other bit A decision for this bit is determined by the sign of the real part of the component. For 
> 0, we use feedback In the sense that the hard decistons for all bits /r* < /rare used to subtract the MAI from the k- 
th component of the output The received signal amplitudes are required for this. Finally, a hard decision for the /r-th 
bit is made. 

ss [0014] SubtmcXwe interference cancellation estimates the transmitted bits in order of decreasing received signal 
powers and hence requires these to be estimated. Initially, the bit estimates provided by a conventional detector are 
employed. For a given bit, the most recent estimates of all other transmission bits are used to generate an estimate 
of the MAI from which it suffers. This interference is subtracted from the signal and an updated estimate of the given 
bit determined by the sign of the real part of the remaining signal. 

40 

OBJECT OF THE INVENTION 

[001 5| The present Invention seeks to provide a simple to Implement base station receiver structure which possesses 
improved symbol detectton characteristics. 

45 

STATEMENT OF THE INVENTION 

[0016] In accordance with a first aspect of the Inventton. there is provided a radto communications system wherein 
an equivalence relation is defined for a set of users. This equivalence relatton is used to group the users into equivalence 

so classes. If CDMA techniques are used, the equivalence relation may be defined according to the temporal infonfnatton 
contained In the partial cross-correlations between the PN codes. If an antenna array is Incorporated into the receiver 
structure at the base statton. the equivalence relation may be defined according to the spatial signatures of the users. 
Alternatively, the equivalence relatton could be defined using both temporal and spatial information. 
If the vectors 25 and are the spatial signatures corresponding to users /and /. a beamforming term for these 2 users 

ss may be defined as (x^/ [jjfjl) (x/ User classes may be determined so that if user / and user y are not in equivalent 
classes, the absolute value of the beamforming term Is small In some sense. Such terms can then be replaced by zero, 
which altows the multiple access Interference between users in different equivalence classes to be ignored. This en- 
ables a reductton In complexity due to the beamforming. 
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[0017] In the case of dispersive multlpath propagation, it may be possible to define more than one spatial signature 
per user at a particular instant. User equivalence classes may be defined as above if one multlpath component (or 
spatial signature) per user is chosen. For example, the first or the strongest multlpath components could be selected. 
The chosen spatial signature for user / will be denoted ago- 
s [0018] Users nnay be grouped into spatial equivalence classes by choosing a beam pattern threshold 8. which lies 
In the following range: O ^ G ^ 1 . User p and user /are said to be close in bearing if: 

llt^'pyiUPVlWf ^« 

10 

[001 9] A first relation V between users is defined by: 

user p . user / if (xj ||ir,.d|f ^ ^ ; 

IS 

and a second relation between users is defined such that: 

user p - user / if there is a finite sequence of users with Indices p^ p^ po., whereby: 

user p s user p^ = user p^ « ... = user Pq = user /. 

[0020] This second relation is an equivalence relation and can be used to determine equivalence classes. 
[0021] Equivalence classes may be used to reduce the complexity of bit detection when using multi-user detection 
techniques. These techniques can be applied within the classes, rather than to all users at once. 

25 

BRIEF DESCRIPTION OF THE FIGURES 

[0022] In order that the present invention can be nriore fully understood and to show how the same may be carried 
into effect, reference shall now be made, by way of example only, to the Figures as shown in the accompanying drawing 
30 sheets wherein: 

Figure 1 is an algebraic representation of a prior art CDMA antenna array system; 

Figure 2 is an algebraic representation of CDMA multi-user detection applied to all users in the CDMA receiver of 
Figure 1; 

3S Figure 3 shows a prior art CDMA antenna array receiver; 

Figure 4 shows a four element unifomi linear antenna array; 

Figure 5 shows a timing diagram for 3 users in a direct sequence signal processing system; 
Figure 6 shows a set of user equivalence classes; 

Figure 7 is an algebraic representation of CDMA multi-user detection applied to user equivalence classes; 
40 Figures 8 a, b, c show a graphical comparison of performance in terms of bit error rate versus number of users for 

various detectors under single path fading propagation conditions at 10 dB signaHo-noise ratio (SNR); 

Figures 9 a, b, c show a graphical comparison of performance In terms of bit error rate versus nunriber of users for 

various detectors under multlpath (J = 3) fading propagation conditbns at 10 dB SNR; ); 

Figures 10 a, b, c show a graphical comparison of performance in terms of bit en-or rate versus number of users 
4S for varbus detectors under additive white Gaussian noise propagation conditions at 10 dB SNR; 

Figures 11 a, b, c show a graphical comparison of performance in terms of bit en-or rate versus number of users 

for various detectors under single path fading propagation conditions at 7 dB SNR; 

Figures 12 a, b, c show a graphical comparison of performance in terms of bit error rate versus number of users 
for varbus detectors under multlpath (J = 3) fading propagation conditions at 7 dB SNR; 

so Figures 1 3 a. b. c show a graphical comparison of performance in terms of bit error rate versus number of users 
for various detectors under Additive White Gaussian Noise propagation conditions at 7 dB SNR; 
Figures 14 a. b, c show a graphical comparison of perfonnance in terms of bit error rate versus number of users 
for varbus detectors under single path fading propagation conditions at 4 dB SNR; and. 
Figures 15 a. b, c show a graphical comparison of performance in terms of bit error rate versus number of users 

55 for various detectors under multlpath (J = 3) fading propagation conditions at 4 dB SNR; and, 

Figures 16 a, b, c show a graphical comparison of perfornnance in terms of bit error rate versus number of users 
for various detectors under Additive White Gaussian Noise propagation conditions at 4 dB SNR. 
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[0023] Note: In Figures 8-1 6, reference (a) refers to a single antenna receiver (array size M=^ ), reference (b) refers 
to a multiple antenna array receiver (M=4) an6 reference (c) refers to a multiple antenna array receiver employing the 
equivalence classes approach. 

5 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0024] There will now be described, by way of example, the best mode contemplated by the inventors for carrying 
out the Invention under multipath and single path conditions. In the foltawing description, numerous specific details are 
set out in order to provide a complete understanding of the present Invention. It will be apparent, however, to those 

10 skilled In the art, that the present invention may be put into practice with variations of the specific. 

[0025] Referring now to Figure 4, a uniform linear antenna array of a telecommunications base station receives 
signals from an outstation (only one is shown for reasons of clarity) within the area of coverage of the base station. 
Figure 5 shows a timing diagram for three users employing direct sequence spread spectrum signals and binary phase 
shift keying. For a number, P. of users (mobiles), the />th mobile will be able to generate a binary data sequence dp 

IS (t), Each symbol has a period of Tg seconds. Each bit for user p is multiplied by the FN code of the user, 

Cp(t), whfeh is a sequence of "chips" with values ±1 , each with period The processing gain, L. is determined by the 
ratfo: L = 7^7^ The transmitted radio frequency signal for the p-th user, s^fl), is given by: 

Sp(t) = Apdp(t)Cp(t)sin(2nft) 

where /denotes the radio frequency of the carrier (with associated wavelength X) and >Ap denotes the signal amplitude. 
Uplink transmissions are taken to be asynchronous, the transmitted signal may be subject to fading and multipath 
propagation. The use of the PN code Cpft) permits the base station to resolve multipath components separated by 7^ 
25 seconds. 

[0026] In the following, the general case of a base station antenna array receiver, operating in a multipath fading 
environment, is considered. 

[0027] By way of example, consider a base station receiver with three uniform linear an-ays. each array consisting 
of Af antennas with spacing A/2 and each covers a single 120* (2?i/3 radians) sector. In this case, the (normalised) 
30 response of the array to a unit impulse from a bearing e is the vector a(6) whose m-th component is equal to: 

exp(l(m ' 1) n sinfi)) 

35 [0028] Using an antenna spacing of X/2 ensures that the mapping: 6 => a(e) Is one-to-one for 9 e n/3, jc/3J radians , 
whilst the scaling factor ensures that 0^(9)5(6)=: 7. 

[0029] Consider a scenario where multipath propagation exists, and each multipath arrives from a single bearing. 
The baseband signal received at the base statbn antenna array can be expressed as: 

40 



45 where Op jp = ApXp jp; Xp ,p is the channel attenuatfon, ^p jp denotes the received signal phase shift; tp^ denotes the 
time delay for the ;p-th multipath component of the p-th user; and, v(t), a complex vector, represents additive noise. 
[0030] To simulate more general Rayleigh fading multipath channels, which permit each multipath to be scattered In 
angle, the model described In [Salz94] (J. Salz and J.H. Winters, "Effect of Fading Correlation on Adaptive Arrays In 
Digital Mobile Radio". IEEE Trans Veh Tech, Vbl 43(4), November 1994, ppl 049-57) is used. This nrxxJel takes Into 

BO account correlations between antenna elements, the antenna spacing and the signal scattering angles. The quantity 
expff WW a(®Pd^/y replaced by: 

55 

Where Cm «s C/V (0,1). a complex normal random variable, 1 Sm^/lf;ii^ 1 ^/n^W,arethe normalised eigenvectors 
of E[Xpjp (t)Xp^ (0^; and. ®^ 1 ^ m ^ M, are the corresponding eigenvalues. Taking these factors into account, the 
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received baseband signal can be expressed as: 



[0031] To detect the yp-th multipath component for the n-th brl of the p-th user, the continuous received vector r^y is 
passed through a temporal filter matched to the PN code CfJ[t) and a spatial filter matched to a scaled version of the 
70 array signature Xp^(t), which can be expressed as: 

15 [0032] This yields a discrete output sampled at the bit rate. Since a single message may contain many bits, the 
message can be divided into groups of N bits, where A/ is a fairly small integer, e.g. /V ^ 1 0, and NP symbols (bits) can 
be decorrelated at a time. For simplicity, it is assumed that Jp=:J,0^p<R\i dpf^^(n) is the n-th transmitted bit from 
the p-th user and y^^^ is the corresponding output detected on the /-th path, 0^r\< N,0^p < P,O^J<J, then: 



20 



25 



30 



40 



45 



SO 



ss 



Y:=(RDATd + z)/L, 

where /? is a positive definite Hermitian matrix of size NPJx NPJ, whose (nPJ+ pJ+l n'PJ+ p'J+f )-th entry is defined 
as; 

Here, x„^j is the discrete time representation of Xp/t); 



P(n.pj)(n',p'4') '® partial cross-correlation for the PN codes of the n-th bit from user p on the /-th path and the n- 
th bit from user p' on the /-th path, assuming P(n,pj)(nj>j) = ^ 
D is a diagonal matrix whose entries are the fading amplitudes Il2&,,pj|; 
>lis a diagonal matrix with diagonal entries A^n,pj)(n4>j) = 
35 TO is a stretched version of d in which each component of d is repeated J times; and, 

z represents noise and E[2^J = o^A, assuming that there is Gaussian noise at the antenna outputs. 

[0033] To combat the effects of fading and to reduce the dimension of the problem from JNP to NP, the outputs are 
multiplied by their corresponding amplitudes and the multipath components are summed to obtain an output 



If 



and 



w=(UL)(AT)^RDAT6 + VL{AT)^Z. 



S = (ATfRDAT\ 



n = OA)(Alfz 



then: 

w=(S/L)d+ri 
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[0034] The matrix S/L represents the transformation from the vector of transmitted bits d to the processed received 
signal w. Vector w will be an approximation to (AT^DATd )Nh\ch can be written as Pd, where P is a diagonal matrix 
whose entries are the sums of the average received powers of muttipath signal components weighted by the fading 

amplitudes. 

5 [0035] The vector d Is estimated using the processed received signal jv. It is possible to estimate d using an alter- 
native signal where: 

w* = (\/L)(AlfDRDATd + {ML) (Alfoz, 

10 

If: 

S = (AjfoRDAT] 

IS 

and: 

n' = 0/L)(ATfDz; 

20 

then: 

(^/L)S'd'hn\ 

25 

[0036] Estimating S' requires less computation than estimating S as no normalisation of array signatures is required. 
Also, S* is Hermitian and positive definite, which is an advantage when applying multi-user detection requiring matrix 
inversion. In this case a Cholesky decomposition can be used to obtain a lower triangular matrix L. where S'^LL^, 
instead of a more general^ decomposition, where^ Is a lower triangular matrix and ^ is an upper triangular matrix. 
30 [0037] However: 

Elrf(n')'^] = (a^/L){S/Ly. 
35 and thus depends upon the fading amplitudes, whereas: 

Eimfj ^(</A)($A): 

40 

where § • = (ATPBAT; and so does not. 

[0038] The fading suffered by jV is more exaggerated than that suffered by w. The bit error rates obtained by esti- 
mating df rom jV are generally much higher than those obtained using w. 

[0039] TTie problem considered in this document arises in the estimation of given the received signal: 

vy= (S/L)d+ n 

[0040] The present approach provides a way of reducing the calculation involved in applying multi-user detection 
so techniques after beamformlng. A single large problem Is split Into a number of smaller problems. It requires a partition 
of the users into equivalence classes. The description given betow Is based on the use of spatial equivalence classes. 
[0041] In the case that N, the number of bits per user, is small; for 0 ^ p < P, 0 ^ < J, it is reasonable to suppose 
the spatial (array) signatures are constant over the time inten^al in question and if it is set such that: 

an equivalence relation on the set of users can be defined based on the spatial signatures ^ where 0 ^ p < P. 
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[0042] Equivalence classes are determined so that, if user p and user / are not equivalent, the absolute value of 
(i"p.</llSi),oll) C5.</lla.oll) is sn^all In some sense. In this case it is reasonable to assume that (Jf^p/ IL^jJI) (Jg/llJ(i.jfl) is 
also'smalf, 0 S / / ^ J. By replacing each such small temi (in fl) by zero, the multiple access Interference between 
equivalence classes can be ignored. This reduces the complexity of multi-user detection applied after beamforming. 
5 The problem is thus split Into a number of small calculations. Figure 6 depicts the partitions of users into equivalence 
classes. 

[0043] The users are sorted into equivalence classes based on the beam patterns according to a beam pattern 
threshold. G, where 0 ^ 9 ^ 1 . User p and user /are close in bearing if: 

[0044] If user pand user / are close in bearing, they are In the same equivalence class. Further, if there exists a finite 
sequence of users, say user p, user qo, user q^, ....user q„ user /, such that any two adjacent users In the sequence 
are close in bearing, user pand user /are in the same equivalence class. Therefore, it is possible that two users in the 
same equivalence class are not close in bearing. However, users in different equivalence classes cannot be close in 
bearing (according to 0). 

[0045] Specifically, a relation V between users is defined by: 

user p « user ||(x"^^/||xJ)(V||xj|fl'^ ^ e. 
[0046] A relation ''^* between users is defined such that: 

user p~ user / if there is a finite sequence of users with indices p^ P|, .... pQ, whereby: 



2S 



3S 



40 



4S 



user p ss user p^ « user p, s= ... = user p^ « user /. 



[0047] Note that the relation » is symmetric and reflexive but not transitive and this is therefore not an equivalence 
30 relation. The relation - Is symmetric, reflexive and transitive and is hence an equivalence relation. By increasing the 
beam pattem threshold 6 we can decrease the number of users in the equivalence classes. If there are a large number 
of users who are fairly equally but closely spaced in bearing, a fairly high threshold should be choosen to obtain a 
practical number of equivalence classes. If 0 = 0 then only one class can be obtained. 
[0048] For a number, E, of equivalence classes and a number of users in the e-th class, then 



can be written as a direct sum, namely: 



wherer* is used to represent the NbWs transmitted by the users in the e-th class. The problem of multi-user detection 
is split into E sma//er problems by considering the restriction Xo<r,^^e^E, which restriction to the e-th equivalence 
so class is denoted by: 

^ [0049] The vectors and matrices contained in this equation model the signal components for the users in the e-th 
equivalence class only. The effects of the other users in the other (E-1 ) classes are nrKxJelled as noise and are included 
In the noise vector n^. 
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[0050] When a signal is transmitted from a base-station to a mobile radb station and vice versa, the signal is reflected 
from buildings and the like so that multiple images of a desired signal arrive at different times. This is known as multipath 
propagation. In addition, transmitted signals are subject to fading. For completeness and simplicity, the situation in 
which no reflections occur (the single path case), and in which Rayleigh fading is absent, will be discussed below. 
5 [0051] For the special case of single path transmission without scattering, the number of multipath components, 
is set to one for all users and a single bearing can be assigned to each user the bearing of the p-th user will be Bp, 
Thus, the received base-band signal at the base station antenna array can for this case be written: 



a special case of equation (1 ). A signal scenario corresponding to equation (3) is simply a special case of equation (1 ) 
^5 and may be processed following the same procedure described above for the case of equation (1). Alternatively it is 
simpler to work with the output: 

y=(Bl)>^flf+z/l 

20 

rather than with jy which Is obtained from ^ by premultiplying by 4". 

[0052] There are a number of alternative approaches to determining equivalence classes for a set of P users. Quan- 
tifying the relationship between users and choosing a threshokJ value for this measurement enables an equivalence 
relation (and hence equivalence classes) to be defined. The relationship between users may be measured in various 
25 ways. Given one choice of measurement, the threshold may be varied. For example, the measurement can based 
purely on the correlation of the vectors xy^ and Xp^, as in equation (2). Some other ways in which the relationship 
between users can be quantified are given bebw: 

(a) The left-hand side (LHS) of equation (2) may be scaled by the cross-correlatfon value of the two PN-codes 
30 corresponding to the 0-th multipath of user / and user p for the A/ symbols that are being demodulated. This cor- 
responds to averaging the quantity (9(n,h0)(nj>,0)^Q over the N symbols. The code correlation levels are likely to 
change more rapidly than the users* spatial signatures, so the equivalence classes may need to be updated more 
rapidly to reflect this. 

(b) The threshold operation in equation (2) may be modified to incorporate the cross-correlation values of all Jf 
35 multipath vectors for user / and Jp multipath vectors for user p. One option is to sum the (J/ Jp) quantities ||(x"p/ 

ll^f^JI) (£yll£,kll)IP possible values of the two users multipath component indices /and k. The result can 

be normalised by the scalar (Jf Jp) and used to replace the LHS of equatton (2). 

(c) An alternative to option (b) is to evaluate the (J{ Jp) correlations and retain only the largest corretatbn value. 
This value may be used to replace the LHS of equation (2). 

40 (d) The measurement may be based purely on the cross-correlation value of the two PN-codes corresponding to 
the 0-th multipath of user / and user pfor the /V symbols that are being demodulated. 

(e) An alternative to option (d) is to evaluate the (Jf-Jf^ PN code correlations and either average over the values 
or retain only the largest correlation value. This value may be used to replace the LHS of equation (2). 

4S [0053] Figure 7 shows a multi-user detectk>n an-angement where data is separated into classes. The components 
of the processed received s]gna\ ws(SA)dfn are partitioned into classes according to user equivalence classes. For 
the e-th user equivalence class, a submatrix S^of the spatial and temporal correlation matrix S is formed. Multi-user 
detectton techniques based on the approximating equatk)n jiv^ s {S^)d^*n^ are then emptoyed to obtain the estimated 
data. 

BO [0054] Each submatrix of Sylekls a good approximation to the effects of signal fading and of pseudo-noise coding 
and decoding on the signals transmitted by the users in its corresponding class. Mutti-iiser detection is thereby restricted 
to these classes, whereby only multiple access terms of little significance are discarded. This results In only a slight 
degradation In the performance of multi-user detectors. 

[0055] The step of combining the multipath components is optional: estimates of d, Adox DATd can be made, using 
55 jvor y If ^ is used, multi-user detection based upon appropriate submatrlces of R/L wouki be applied to the classes. 
If the equivalence classes are selected appropriately, for each class the appropriate submatrix of R/L would yield a 
good approximation to the temporal and spatial cross-correlations between the signals transmitted by the users in the 
class, up to fading effects. 
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[0056] Simulations were carried out in order to compare the bit error rates achieved by three approaches to bit 
detection. The first approach employed a single antenna and standard bit detection techniques; the second approach 
employed beamf orming by means of a uniform linear array of Mantennas, M> 1 , followed by the applicatbn of standard 
bit detection techniques; and the third approach employed beamforming as for the second approach, partitioning the 
5 users into equivalence classes based on their spatial signatures, followed by the application of standard bit detection 
techniques, restricted to the classes. 

[0057] The simulations were based upon a scenario of up to 33 users per cell, with an array of four antennas for the 
beamforming and Gold codes of length 31 were used. The number of bits per user considered at one time was set at 
four (N = 4) and the beamforming threshold for third approach was set at 9 = 0.7. Noise samples at the antenna outputs 

70 were assumed to be Gaussian, and signal to noise ratios of 4, 7 and 1 0dB were selected. Five methods of bit detection 
were simulated. The simulations were based on three channel types: 1 ) additive while Gaussian noise (AWGN) chan- 
nels; 2) single path (J= 1) Rayleigh fading channels, and; 3) multipath (J= 3) Rayleigh fading channels. Graphical 
representations of these measurements are displayed as bit error rates versus number of users in Figure 8-16. 
[0058] To model the AWGN channels, single path transmission was assumed. The time delays (measured in chips) 

IS were uniformly distributed over [0,31]. The channel amplitudes and attenuations were all set to one and the phase 
shifts uniformly distributed over [0. 2ic]. The bearings of the users were uniformly distributed over [-it/3,n/3]. 
[0059] To model the Rayleigh fading channels the time delays of first multipath components were uniformly distributed 
over [0.31]. The central bearings of first multipath components were uniformly distributed over [-ic/3,7i/3). For each user, 
central bearings of other multipath components were uniformly distributed about the first with a spread of one radian. 

20 It was assumed that each channel tap arrives from Its central bearing with uniform angular power density over a spread 
of 0.33 radians. For the single path Rayleigh fading channels the average received power per user was set to one. For 
multipath transmission; for each user, relative path delays were set to one chip; and average received powers per user 
to 1,0.5, and 0.25. 

[0060] Figures 8a-1 6a (depicting the first approach) confirm that multi-user detection yields lower bit error rates than 
25 conventional detection and hence Increases capacity. This increase in capacity becomes more significant as the signal- 
to-noise ratio increases or as fading effects decrease. 

[0061] For the first approach, in general, subtractive interference cancellation yields the lowest bit error rates. For 
very high numbers of users, the decision feedback decorrelator may perform best. Subtractive interference cancellation 
provides almost single-user performance in the presence of severe fading. 
30 [0062] Reducing multiple access Interference by beamforming enables all detectors to achieve lower bit error rates 
and hence Increases capacity (the second approach, see Figures 8b-16b). Significant capacity gains can be achieved 
by using beamforming followed by multi-user detection if fading Is not too severe and signal-to-noise ratios are high 
enough. 

[0063] The bit error rates achieved by multi-user detectors become closer as multiple access Interference is reduced. 

OS Subtractive Interference cancellation provides the lowest bit error rates, giving almost single-user performance. 

[0064] Beamforming followed by multi-user detection applied to spatial equivalence classes (the third approach, see 
Figures 8c-16c) yields slightly higher bit error rates than beamforming followed by multi-user detection applied to all 
users at once. As the signal-to-noise ratio decreases or as fading effects Increase, the difference in bit error rate given 
by the two approaches becomes less significant. This approach Is potentially much less computationally complex than 

<o the second approach. 

[0065] For the third approach, subtractive interference cancellation provides the lowest bit error rates and tends 
towards single-user performance as signal-to-noise ratio decreases or as fading effects increase. " " 

[0066] It has been shown, for the reverse link of a cellular direct sequence code division multiple access system, 
that by using mutti-user detection (Moshavi96] or antenna array receivers [Naguib94]. bit error rates considerably lower 

45 than those provided by the conventional detector can be achieved. The results depk:ted in Figures 8a.b-16a,b (first 
and second approaches) confirm these findings for both single path AWGN channels and Rayleigh fading single and 
multipath channels and show that Initial beamforming followed by multi-user detection further decreases bit error rate 
and hence Increases capacity, but at a cost in complexity. The results depicted in Figures 8c-16c (the third approach) 
indicate that initial beamforming followed by multi-user detection applied to spatial equivalence classes is worthy of 

so consideration. 

[0067] In general, the conventional detector gives the highest bit error rates and subtractive interference cancellation 
gives the lowest bit error rates. The decision feedback detector may perfomi best for large numbers of users if beam- 
forming is not used The MMSE and decision feedback detectors perfomi better than the decorrelator. 
[0068] Significant capacity gains can be achieved by using beamforming followed by multi-user detection if fading 
ss Is not too severe and signal-to-noise ratios are high enough. The second approach gives the highest capacity, but is 
the most costly to implement If the number of antennas at the receiver is Increased, beamforming techniques can be 
refined and the sizes of the equivalence classes can be reduced. By restricting to equivalence classes, the computa- 
tional burden can be reduced. By choosing the beamfomiing threshold wisely, only multiple access Interference terms 
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of little significance are discarded and hence only a slight degradation In the performance of multi-user detectors is 
observed. 

[0069] When restricting to equivalence classes, subtractlve Interference cancellation generally provides the best bit 
error rates. It is more computationally costly than conventional detection; but much less costly than using the decor- 
s relator, the MMSE detector implemented in this study or the decision feedback decorrelator, all of which are based on 
matrix inversion. However, like the matrix inversion implementatfon of the MMSE detector and the decision feedback 
decorrelator, it requires a good deal of informatkDn, namely, the timing and codes of all the users, and the relative 
amplitudes of the incoming signals. 

[0070] The computational complexity incurred by beamforming depends upon the algorithm selected for determining 
10 the spatial filters [Thompson96] (J.S Thompson, P.M. Grant and B. Mulgrew. "Performance of Antenna Array Receivers 
for CDMA", Pfoc. IEEE Globecom Conf. London (UK), November 1996, pp570-4.). Different algorithms vary in their 
complexities. However, multi-user detection applied to spatial equivalence classes is potentially far less complex than 
multi-user detection applied to all users at once. Hence, if fading is not too severe and signal-to-noise ratios are high 
enough, the present invention can provide, using one of the simpler forms of multi-user detectfon, the best compromise 
15 between performance and computational cost. 



Claims 

20 1. A radio communications system wherein an equivalence relatton is defined for a set of users based on a set of 
beam-forming terms to group the users according to their spatial signatures. 

2. A system according to claim 1 wherein the beam-forming terms are aW(e^a(e^. 

25 3. A system in accordance with claim 2 wherein the user classes are determined so that, if user / and user y are not 
in equivalent classes, the absolute value of a^(0j)a(8^ is small in some sense. 

4. A system in accordance with claim 2 wherein the absolute value of a^(6/)a(ej) is then replaced by zero, which 
allows the multiple access Interference to be ignored between equivalence classes. 

30 

5. A system in accordance with claim 2 wherein the users are sorted into equivalence classes based on the beam 
pattems, where the beam pattern has a threshoki 9 which lies in the following range: O < 0 ^ 1 ; in accordance 
with the folbwing relations: 

35 a first relatksn V between users, which is symmetric, reflexive and non-transitive, and is defined by: 

user / « user ; if ||a"(ej)a{ep||^^ 0; 

40 and 

a second relation between users, which is symmetric, reflexive and transitive, and is defined by: 

user / — user I 

45 

provided that there exists a finite sequence of users with indices it -^ht ^"^^ 

user / s= user » user = ... » user Iq user / 

so 

6. A system according to claim 1 , wherevi the beam-forming Xems are {jf^pJlix^JUQ^f^f/liXfJ]), 

7. A system in accordance with claim 6, wherein the user classes are determined so that, if userpand user /are not 
in equivalent classes, the absolute value of {J^pj J) {Xf^ is small in some sense. 

55 

8. A system in accordance with claim 6, wherein equivalence classes can be determined in this situation if it can be 
assumed that if the quantity {j^^J ||Xp J) {x^J Iji/J) is small, the quantities (xHp/ lU^^^jl) Osy/ where 0 ^ y, 
/ ^ J, are also small. This permits all the multiple access interference between two users in different equivalence 
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IS 



classes to be ignored. 

9. A system in accordance with claim 6. wherein user p and user /are sorted into equivalence classes based on the 
beam patterns according to a beam pattern threshold, 9, which lies in the following range: O s 0 ^ 1 ; and, 
wherein user p and user / are close in bearing if: 



10 10. A system in accordance with claim 6, wherein the users are sorted into equivalence classes based on the beam 
patterns, where the beam pattern has a threshold G which lies in the folbwing range: O < G ^ 1; In accordance 
with the following relations: 



a first relation V between users is defined by: 

user p» user /« l(i"M\)(sM\f ^ 

and 

20 a second relation '~' between users is defined such that: 

user p user / if there is a finite sequence of users with indices p^ p^, Pq, whereby: 

user p =5 user p^ s user p^ « ... « user p^ « user /. 

11. A system according to claim 1 , wherein the modulation scheme Is a code division multiple access scheme. 

12. A method of operating a radio communications system wherein an equivalence relation Is defined for a set of users 
30 based on a set of beam-forming terms wherein the users are grouped according to their spatial signatures. 

13. A method according to claim 12 wherein the beam-forming terms are a''(e^a(ep. 

1 4. A method according to claim 1 2. wherein the beam-forming terms are Q^pJ \\XpJ\) QCfJ IteJU 

35 

15. A method according to claim 12, wherein the modulation scheme is a code division multiple access scheme. 



40 



45 



so 



55 



12 



EP 0 926 913 A2 




^ C 

.S Ml 
<0 

— ** 

2 X 
o c O 

I- c 

(0 is 
■5 I 

CO 



o 


NRJ 


_J 


A 

• 




z 



I— 

< 



CO 

<^ CO 

OQ c 



i8 



C CO 

♦? o 



.2 o 

3 « S 

so $ 

CQ 



gco 



<T3I 


• 


All Users 

NP->NP 


^1 





.= o. e 



P 

w CD 
S 



Z 



x»| 



13 



EP 0926 913 A2 




14 



EP 0926 913 A2 




UserO 
Userl 
User 2 



BitO 



Bit 3 



Bite 



Bit 9 



Bit1 


Bit 4 


Bit? 


Bit 10 



Bit 2 


Bits 


Bits 


Bit 11 



r4—4 — H 



0 M2TS tj^JTs 

Fig. 5 

Equivalence Classes of Users 

I 



Receiver Array | 




Fig. 6 



16 



EP 0 926 913 A2 



Partition Components Multiuser Detection Bit 

According to Based on Estimates 

Equivalence Classes Submatrices of S/L 



1st Class 


^1 


1st Class 


3, 


NP •> NP, 




NP,->NP, 




• 


• 


• 


• 


e-th Class 


We 


e-th Class 


A 


NP->NP, 




NFJ->NP^ 




• 








E-th Class 




E-th Class 




NP->NfJ 




NPg->NPg 





Fig. 7 

16 



EP 0 926 913 A2 




EP 0926 913 A2 




EP 0 926 913 A2 




EP 0926 913 A2 




20 



EP 0 926 913 A2 




EP0 926 913 A2 




EP0 926 913 A2 




EP 0 926 913 A2 




EP 0 926 913 A2 




EP 0 926 913 A2 



I 



E ^ o 
c 



ill 



n I 



till! 



O 

o 



e^oy JQJLJ3 )!g 



26 



EP 0 826 913 A2 



1 • I s 

t p ? + 



OTD t ^ O 

c 



Hi 



' * ■ ■ ■ 



o 

a^oy JQJJ3 ^la 



o 
ci 



27 



EP 0 926 913 A2 




EP 0 926 913 A2 




EP 0926 913 A2 




EP 0 926 913 A2 




EP 0 926 913 A2 




32 



EP 0 926 913 A2 




JCUJ3 )!Q 



33 



EP0 926 913 A2 



2d 




1 1 ■ ■ ■ ■ ■ 



-L 



i I ■ I i 



ID 



o 

CO 



to 



O (U 



z 



IT) 



o 
o 



o 



epa J0JJ3 ^18 



o 
o 
o 
d 



34 



EP 0 926 913 A2 



t a M 



JD 



1 1?! 

i 'i i 

I * • 



1 is I 



'•1! 



I I I I I ■ 



in 



CM 



1/1 

(VI 



z 



13 



o 

5 
o 



BjDy JQJJ3 «g 



35 



EP0g26 913 A2 




O 

e^oy JQJJ3 m 



o 
o 



36 



EP 0 926 913 A2 




9)Dy JQJJ3 l!9 



38 



EP 0 926 913 A2 




O 



o 
o 



39 



EP 0 926 913 A2 



I 

c 

JQ 

in 



fM 



3 



in 



I I I i I 

r- 
O 

a^oy jajJ3 ^ig 



5 



40 



EP 0 926 913 A2 



in 




apy J0JJ3 ^ig 



41 



EP 0 926 913 A2 




EP 0926 913 A2 




(19) 



J 



EuropSlsches Patentamt 
European Patent Office 
Office europeen des brevets 



liiiliilliilllili 

(11) EP 0 926 913 A3 



(12) 



(88) Date of publication A3: 

21.06.2000 Bulletin 2000/25 

(43) Date of publication A2: 

30.06.1999 Bulletin 1999/26 

(21) Application number: 98203975.2 

(22) Date of filing: 25.11.1998 



EUROPEAN PATENT APPLICATION 

(51) Intci7: H04Q 7/36. H04B 7/26 



(84) Designated Contracting States: 

AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 
MC NL PT SE 

Designated Extension States: 
ALLTLVMK ROSI 

(30) Priority: 24.12.1997 GB 9727352 

(71) Applicant: NORTEL NETWORKS CORPORATION 
Montreal, Quebec H2Y 3Y4 (CA) 

(72) Inventors: 

• Hassell Sweatman, Catherine Zoo Wollaston, Dr. 
Mayfleld Road, Edinburgh, EH9 3JL (GB) 



• Thompson, John Scott, Dr., 
University of Edinburgh 

Mayfleld Road, Edinburgh, EH9 3JL (GB) 

• Mulgrew, Bernard, Dr., University of Edinburgh 
Mayfleld Road, Edinburgh, EH9 3JL (GB) 

(74) Representative: Humphrey-Evans, Edward John 
Intellectual Property Law group, 
Nortel Networks, 
London Road 

Harlow, Essex CM17 9NA (GB) 



(54) Cellular radio communication system with user partioning by antenna beam forming 



(57) The present Invention relates to radio commu- 
nication antennas and In particular relates to multi-user 
detection for CDMA antenna array receivers. A reverse 
link (mobile to base) of a cellular direct sequence code 
division multiple access system is provided which in- 
creases system capacity whilst avoiding system com- 
plexity. The invention provides a method of signal de- 
tection which is based on the partitioning of users Into 



equivalence classes by beam-forming with subsequent 
standard multi-user detection techniques applied within 
each class. Computation techniques are simplified and 
processing time can be considerably reduced Further- 
more, the capacity of a conventional receiver may be 
improved by employing multi-user detection techniques. 
These techniques are particularly useful for systems 
which Involve large numbers of users. A method of op- 
eration is also disclosed. 
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